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ABSTRACT 

We determine the contribution of stars in galaxies, intracluster stars, and the intracluster medium 
to the total baryon budget in nearby galaxy clusters and groups. We find that the baryon mass 
fraction (/t = f2f,/0m) within r5oo is constant for systems with iWsoo between 6 x 10^^ Mq and 
1 X 10^^ M0. Although fb is lower than the WMAP value, the shortfall is on the order of both the 
observational systematic uncertainties and the depletion of baryons within rsoo that is predicted by 
simulations. The data therefore provide no compelling evidence for undetected baryonic components, 
particularly any that would be expected to vary in importance with cluster mass. A unique feature 
of the current analysis is direct inclusion of the contribution of intracluster light (ICL) in the baryon 
budget. With the addition of the ICL to the stellar mass in galaxies, the increase in X-ray gas mass 
fraction with increasing total mass is entirely accounted for by a decrease in the total stellar mass 
fraction, supporting the argument that the behavior of both the stellar and X-ray gas components is 
dominated by a decrease in star formation efficiency in more massive environments. Within just the 
stellar component, the fraction of the total stellar luminosity in the central, giant brightest cluster 
galaxy (BCG) and ICL (hereafter the BCG-flCL component) decreases as velocity dispersion (a) 
increases for systems with 145 km s^^ < cr < 1026 km s~^, suggesting that the BCG-I-ICL component, 
and in particular the dominant ICL component, grows less efficiently in higher mass environments. 
The degree to which this behavior arises from our sample selection, which favored systems with central, 
giant elliptical galaxies, remains unclear. A more robust result is the identification of low mass groups 
with large BCG-I-ICL components, demonstrating that the creation of "intracluster" stars does not 
require a massive cluster environment. Within rsoo and r200j the BCG-I-ICL contributes on average 
40% and 33% of the total stellar light, respectively, for the clusters and groups in our sample. Because 
these fractions are functions of both enclosed radius and system mass, care should be exercised when 
comparing these values with other studies and simulations. 

Subject headings: galaxies: clusters: general — galaxies:cD, formation, evolution, fundamental pa- 
rameters 



1. INTRODUCTION 

An important step in understanding any astrophysical 
system is a full accounting of its constituent components. 
Over scales ranging from groups to clusters of galaxies, 
measuring how baryons are distributed among their var- 
ious gaseous and stellar phases can provide clues to the 
roles played by AGN feedback, starburst winds, galaxy 
mergers, and tidal stripping (Merritt 1984; Roussel et al. 
2000; Valdarnini 2003; De Lucia et al. 2004; Kravtsov 
et al. 2005; Cheng et al. 2005; Ettori et al. 2006; Monaco 
et al. 2006) in driving galaxy evolution in these com- 
plex environments. Determining the total contribution 
of these baryonic components to the cluster mass is simi- 
larly valuable. While the baryon fraction in clusters may 
be expected to reflect the universal value (see White et al. 
1993; Ettori et al. 2006), baryon accountings in clusters 
to date have fallen short of this expectation (Ettori 2003; 
McCarthy et al. 2006) and uncovered possible intrigu- 
ing trends with cluster mass (Lin et al. 2003). With 
the precise measurement of the universal baryon fraction 
from WMAP (/& = n^/n^ = QAIQ^q^^iI] Spergel et al. 
2006), this shortfall has gained in physical significance. 
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leading to suggestions that physical processes lower the 
cluster baryon fraction relative to the Universe (see, for 
example. He et al. 2005), that there may be significant, 
undetected baryon components (see Ettori 2003; Lin & 
Mohr 2004), or that WMAP underestimates (Mc- 
Carthy et al. 2006). Any investigation along these lines 
first requires measurement of all the baryon components, 
including the hot. X-ray emitting gas and the stars in and 
out of galaxies. 

One baryonic component that past studies were un- 
able to include directly is the intracluster stars (here- 
after the intracluster light or ICL). In recent years, mul- 
tiple studies have detected intracluster light via stacking 
of large cluster samples (Zibetti et al. 2005), deep ob- 
servations of individual clusters (Feldmeier et al. 2004b; 
Gonzalez et al. 2005; Krick et al. 2006), and even identifi- 
cation of individual intracluster stars in very nearby clus- 
ters (Durrell et al. 2002; Feldmeier et al. 2004a; Aguerri 
et al. 2005; CiarduUo et al. 2005; Gerhard et al. 2005). 
The general consensus is now that all clusters have an 
ICL component. The natural extension of these studies 
is placing the ICL in the larger context of cluster and 
group evolution, examining directly whether this com- 
ponent is important in the chemical enrichment of the 
intragalactic medium (Lin & Mohr 2004; Zaritsky et al. 
2004), the total stellar budget, or even in the total baryon 
budget. While recent work suggests that the ICL does 
not contribute much to the overall baryon budget of the 
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most massive clusters (the ones most often observed), 
it could become important in lower mass systems where 
the plasma is less dominant and the apparent discrep- 
ancy between observations and the WMAP baryon frac- 
tion is most severe (e.g. Lin et al. 2003). This issue is 
mostly unexplored, as previous ICL measurements have 
come either from a few systems with little dynamic range 
(Feldmeier et al. 2004a; Krick et al. 2006) or from stacked 
systems for which mass estimators have been crude (Zi- 
betti et al. 2005). 

In our program we include the ICL's contribution to 
the total baryon mass budget directly, explore the de- 
pendence of total baryon fraction on cluster mass, and 
consider the relative importance of the stars in the ICL 
to those in galaxies and of the total stellar mass to 
the total gas mass in the intracluster medium (ICM). 
We focus on a sample of 23 nearby galaxy clusters and 
groups with ICL measurements from Gonzalez et al. 
(2005) that span a wide range in velocity dispersion, 
145 km s""' < o" < 1026 km s"-'^. We determine the 
cluster radii and masses at an overdensity of 500 times 
critical density (rsoo and !M'5oo;§3.1), revisit our previ- 
ous ICL measurements (§3.2), measure the total lumi- 
nosity in the cluster galaxies (§3.3), determine stellar 
masses (§3.4), and obtain gas masses from the litera- 
ture (§3.5). We then derive relationships for the stellar 
and total baryon mass fraction with cluster mass (§4.1), 
examine the detailed behavior of the ICL relative to the 
galaxy luminosity (§4.2), and speculate about the origin 
of the observed trends (§4.3). Throughout this paper we 
assume a cosmology with = 0.27, = 0.73, and 
Ho = 70 km s-\ 

2. SAMPLE DEFINITION AND DATA 

In Paper I (Gonzalez et al. 2005), we present a sample 
of 24 clusters and groups for which we obtained drift 
scan imaging in Gunn i using the Great Circle Camera 
(Zaritsky et al. 1996) on the Las Campanas Im Swope 
telescope. The sample is comprised of systems at 0.03 < 
2 < 0.13 that contain a dominant brightest cluster galaxy 
(BCG) with a major axis position that lies within 45° 
of the direction of the drift scan, east- west, and span 
a range of velocity dispersions and Bautz-Morgan types 
(Bautz & Morgan 1970). We refer the reader to Paper I 
for further details concerning the data and photometric 
reductions. 

In Paper II (Zaritsky et al. 2006), we present velocity 
dispersions for 23 of these clusters and groups. Twelve 
of these a's are based on our own observations with the 
Hydra spectrograph at the CTIO Blanco 4m telescope, 
while the rest are derived from galaxy redshifts in the lit- 
erature. For all but four systems, the velocity dispersions 
are calculated using > 20 confirmed members within 1.5 
Mpc of the BCG. The system with the smallest number 
of confirmed members has 8 galaxies and the lowest a 
(144^57 km s"-'^). Details of the spectroscopic reduc- 
tions and velocity dispersion measurements are given in 
Paper II. The information from both Papers I and II rel- 
evant to the current analysis is presented in Table I . The 
magnitudes in this Table do not include extinction or k- 
corrections, which are listed in a separate column. When 
converting to stellar masses in §3.4 we do include these 
corrections. 



3. ANALYSIS 

To determine the ratio of baryons to total mass for 

each group and cluster in our sample, we 1) estimate the 
total masses and cluster radii, 2) calculate the combined 
luminosity of the brightest cluster galaxy and intracluster 
light, 3) determine the total luminosity associated with 
the other cluster members, 4) convert the luminosities 
of these components into stellar masses, and 5) obtain 
X-ray gas mass fractions from the literature. In §4, we 
discuss the relationships between total system mass and 
stellar mass. X-ray gas mass, and total baryon mass. 

3.1. Cluster Radii and Masses 

To estimate the total baryon fraction from the stellar 
components and ICM, we must first choose a fixed radius 
within which the baryons can be summed and then de- 
termine the baryonic and total masses within this radius. 
We choose r5oo, the radius at which the cluster mass den- 
sity exceeds the critical value by a factor of 500, and the 
largest radius for which the current X-ray data require 
no model extrapolation (Vikhlinin et al. 2006). Mea- 
surements of the ICL typically do not extend even out 
to rsoo, but the light is sufficiently concentrated within 
a few hundred kpc that the model extrapolation to rsoo 
introduces a minimal uncertainty. 

The systems in our ICL sample have measured veloc- 
ity dispersions, but generally lack X-ray data. To obtain 
rsoo and CM500 values for the clusters in our ICL sample 
that are directly comparable to those determined in X- 
ray studies, we therefore derive calibrations of the a- rsoo 
and (T-9Vfsoo relations using the clusters from Vikhlinin 
et al. (2006) that also have published velocity disper- 
sions (Girardi et al. 1998; Wu et al. 1999). The fit to 
the cT-rsoo relation has a slope consistent with the theo- 
retical expectation (1.07 ±0.12 vs 1) and corresponds to 
a cluster with cr = 1000 km having an rsoo of 1.41 
Mpc. The best fit to the cr-fWsoo relationship using the 
Vikhlinin ct al. (2006) sample has a slope consistent with 
the theoretical expectation (3.37 ± 0.53 vs. 3) and cor- 
responds to a cluster with a = 1000 km s~^ having an 
enclosed mass, iAfsoo, of 7.6 x 10^'* M0. A similar analysis 
can be performed for r2oo using the sample of Arnaud 
et al. (2005), but with the limitation that only five clus- 
ters in this sample have published velocity dispersions. 
Fixing the slope to the value derived for the Vikhlinin 
ct al. (2006) sample, the resulting normalization yields 
an r2oo of 2.2 Mpc for a a = 1000 km s~^. We utilize 
this r2oo relation in §4.2.3 when quantifying the radial 
dependence of the ICL. 

For comparison, we also derive rsoo and r2oo using the 
approach of Hansen et al. (2005), directly calculating 
the radius within which the number density of cluster 
galaxies implies that the mass density exceeds the criti- 
cal value by a factor of 200 or 500. If the bias between 
the galaxy distribution and dark matter is not large, this 
approach is expected to yield a reasonable estimate for 
these radii. As in Hansen et al. (2005), we use the Blan- 
ton et al. (2003) SDSS i— band huiiinosity function in- 
tegrated over redshift to compute the field density (see 
Hansen et al. 2005). When computing r2oo and rsoO) we 
impose a faint magnitude limit of mj = 18; varying this 
limit by ±0.5 mag alters the derived r2oo values by ±5%. 
We find that two approaches yield consistent values ex- 
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TABLE 1 
Cluster Sample 
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subtraction, the background subtraction for the galaxy luminosity, and the uncertainty in the faint end slope of the cluster galaxy 
luminosity function. 

^For both the BCG+ICL and galaxies, the quoted absolute magnitudes do not include the corrections for extinction or k-dimming 
listed in column 8. 
'^The estimated uncertainty for the galaxies is 0.15 mag. 

^ The extinction corrections are based upon the Schlegel et al. (1998) dust maps and the k-corrections assume a passive stellar 
population. 

Abell 2405 is a superposition of two relatively poor systems at 11,000 and 27,000 km . While this superposition hinders 
a robust measurement for the luminosity of the cluster galaxies, the BCG+ICL fraction is > 80% within r5oo even if no galaxy 
background subtraction is performed. 

*No published Bautz-Morgan type. Abell 2405 and Abell 3693 each have two discrete redshift peaks. Published Bautz-Morgan 
types do exist for these clusters, but not for the redshift peaks that we are studying. 



cept for clusters with low velocity dispersions (< 400 
km s~^), where the Hansen et al. (2005) approach yields 
larger radii (by as much as a factor of two) . For all anal- 
yses below we use the X-ray derived values; no results 
in the paper change qualitatively if we instead use the 
values from the Hansen et al. (2005) approach. 

3.2. Total Luminosity of BCG+ICL 

A key question in previous discussions of these data 
(Zaritsky et al. 2004; Gonzalez et al. 2005) was whether 
the BCG and ICL are distinct entities. A central result 
of Paper I is that the stellar surface brightness, position 
angle, and ellipticity profiles of these clusters are best 
reproduced when the BCG and ICL are fit as separate 
components. In this analysis we choose to combine the 
luminosity of those two components and refer to the com- 
bined entity as simply the BCG+ICL. There are several 
strong motivations for this choice in the current context. 
Observationally, the separation of the two components 
is difHcult and the combined luminosity is therefore a 
more robust quantity than the luminosity of either sep- 



arately. Using the combined quantity also avoids choos- 
ing among the various ICL definitions in the literature 
(e.g., Feldmeier et al. 2004b; Zibetti et al. 2005), enabling 
a more straightforward comparison between our results 
and other observational studies and simulations. Lastly, 
because we are interested here in a full baryon account- 
ing, the question of whether they reside in the central 
galaxy or in intracluster space is not as critical as in 
other contexts. In later discussions of the ICL relative to 
the galaxy luminosities (§4.2), where the distinction be- 
tween BCG and ICL is relevant, wc reiterate from Paper 
I that the bulk (> 80%) of the luminosity in the cluster's 
BCG+ICL is contained in the ICL. 

We determine the luminosity of the BCG+ICL compo- 
nent as follows. In Paper I we used the GALFIT package 
(Peng et al. 2002) to fit the surface brightness distribu- 
tion in each cluster out to a radius of 300 kpc from the 
BCG. The BCG and ICL were best fit with separate r^/^ 
profiles. Here we use those best-fit profiles to construct 
a 2-D model image from which we determine the fiux 
within circular apertures corresponding to rsoo (in our 
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analysis of the baryon mass fraction; §4.1) and to frac- 
tions of r2oo (in our determination of the BCG+ICL's 
contribution to the total stellar light; §4.2). We esti- 
mate the systematic errors in the measured fluxes us- 
ing the ±lcr systematic uncertainties in the structural 
parameters of the best-fit profiles (these exceed the sta- 
tistical errors due to uncertainties in the sky level; see 
Paper I). The advantages of using the models rather 
than placing apertures on the original data are that we 
minimize contamination from cluster galaxies, interpo- 
late directly over masked regions, and provide a robust 
means of quantifying the systematic uncertainties. We 
use circular rather than elliptical apertures to enable di- 
rect comparison with other studies and simulations. 

An implicit assumption in this work is that our mod- 
els from Paper I, which are determined for r < 300 kpc, 
adequately describe the ICL at all radii. In our models, 
80% of the total fight of the BCG+ICL is contained in 
the central 300 kpc, making it unlikely that we are signif- 
icantly overestimating the luminosity of this component 
with our extrapolation. Moreover, the BCG-I-ICL domi- 
nate the total cluster luminosity within this radius, and 
for our data the potential bias from inclusion of faint, 
unresolved cluster galaxies in the ICL is at most a few 
percent (Gonzalez et al. 2000). Although the uniformity 
of the profiles within 300 kpc argues that extrapolating 
is reasonable, we will underestimate the ICL luminos- 
ity if there is any significant contribution at larger radii 
that is not included in our model. We expect that the 
best method of testing this assumption is via extension 
of current planetary nebulae studies in nearby clusters 
like Virgo (Feldmeier et al. 2004a) to large fractions of 
r2oo- Because of the nature of this study, we are most 
concerned with overestimating the ICL, and large over- 
estimates (> 30 %) are not possible. 

3.3. Total Luminosity of Cluster Galaxies 

We compute the total luminosity of non-BCG cluster 
galax;ies within the same apertures (rsoo and fractions 
of r2oo) used for the BCG-I-ICL. Hereafter, we refer to 
the non-BCG cluster galaxies as simply "cluster galax- 
ies" . After excluding stars using the SExtractor (Bertin 
& Arnouts 1996) stellarity index, we sum the flux of all 
galax;ies fainter than the BCG and brighter than m/ = 18 
lying within the aperture. We then perform a statistical 
background subtraction, and include a completeness cor- 
rection to account for the contribution of fainter cluster 
members (see below) . In cases where the aperture size is 
larger than the width of the drift scan images, we also ap- 
ply a correction to account for the lost area. Specifically, 
we weight the luminosity contribution of each galaxy by 
the inverse of the fractional area lost at that radius (e.g., 
a galaxy would have a weighting factor of two if half the 
area is lost). 

Because the determination of the flux density from 
background galaxies is critical in this type of measure- 
ment, we estimate the background using two indepen- 
dent methods. First, we compute the background galaxy 
density with an annular region located 30' to 60' from 
the BCG. This method, which was also used in Paper 
I, has the advantage that the photometry for the back- 
ground and cluster galaxies is performed in the samc^ 
fashion from the same images. However, large scale 
structure may bias the observed background level even 



at these large radii. To ascertain whether we are suffer- 
ing from such a bias, we also compute the background 
level using the Sloan Digital Sky Survey luminosity func- 
tion from Blanton et al. (2003), integrating over redshift 
and applying mean evolutionary and k-corrections de- 
rived from Bruzual & Chariot (2003) models with the 
Padova isochrones (Bertelli et al. 1994). A comparison 
of the results from the two techniques shows that for 
the ensemble of systems the total cluster galaxy lumi- 
nosities after background subtraction are consistent to 
within 10%. In the subsequent discussion, we apply the 
first method of direct background subtraction because 
this approach relies on fewer assumptions. 

After background subtraction, we include a complete- 
ness correction to account for the luminosity contribution 
from cluster members fainter than the limiting magni- 
tude. For this correction, we adopt the cluster luminosity 
function of Christlein & Zabludoff (2003), a = -1.21 and 
M*,i{ = -21.14, using R- I = 0.82 to convert from R to 
/. Although the value of M* is identical to the i— band 
value for the field from Blanton et al. (2003), the clus- 
ter luminosity function has a slightly steeper faint end 
slope. Using a different value of a systematically shifts 
the total luminosity due to cluster galaxies and thus the 
total stellar light of the cluster. For a = —1, the com- 
puted total luminosity of the cluster galaxies would on 
average be 12% lower, increasing the relative importance 
of the BCG+ICL, but decreasing the stellar contribution 
to the overall baryon budget. The total luminosities of 
the BCG-I-ICL and of the cluster galaxy components are 
presented in Table 1. 

3.4. Total Stellar Masses 

The conversion from luminosity to stellar mass for 
the cluster galaxies and the BCG-I-ICL potentially in- 
curs the most significant uncertainty in the total baryon 
mass biidget. A straightforward theoretical calcula- 
tion, based on PEGASE2.0 (Fioc & Rocca-Volmerange 
1997) for a Salpeter IMF and a 10 Gyr old population, 
yields M/Lj = 4.7, but recent observations indicate that 
real systems exhibit lower M/L. To estimate a mean 
M/L, we use the SAURON results from Cappellari et al. 
(2006), which provide an empirical determination of the 
central stellar M/L using Schwarschild dynamical mod- 
eling of two-dimensional kinematic data. Equation 9 
in Cappellari et al. (2006) provides the luminosity de- 
pendence of the mass-to-light ratio in the /—band. We 
compute a luminosity-weighted M/L for L > 0.25L* 
(the range over which the SAURON relation is deter- 
mined) using the same Christlein & Zabludoff (2003) lu- 
minosity function as in §3.3. We adopt the resultant 
value, M/L = 3.6, for both the cluster galaxies and the 
BCG-j-ICL in our baryon calculation, but note that the 
uncertainty is at least 10% and that M/L varies with 
both mass (Cappellari et al. 2006; Zaritsky et al. 2006) 
and stellar population. An improved accounting for the 
ICL component would require knowledge of the source of 
the ICL stars if the variations in M / L are driven by popu- 
lation, rather than dynamical, differences. An improved 
accounting for the cluster galaxy population would re- 
quire determining the fraction of spiral galaxies in these 
clusters and using a more appropriate M/L for the spi- 
rals. 
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Fig. 1. — Determination of cluster and group baryon fractions within rsoo as a function of sw^soo (bottom axis) and velocity dispersion 
(top axis). Plotted are X-ray gas mass fractions from Vikhlinin et al. (2006) (filled circles) and Gastaldello et al. (2006) (filled triangles), 
and the stellar mass fractions (BCG+ICL+galaxies, open squares) for systems in our sample with masses that overlap the range shown 
for the X-ray studies. All measurements are within rsoo- Overplotted are the best fit relations for the Vikhlinin et al. (2006) sample gas 
mass fractions and for our stellar mass fractions. The WMAP la confidence region for the universal baryon fraction from Spergcl et al. 
(2006) is shown for comparison, and the right-hand axis shows Y, the ratio of the baryon fraction for each component to the universal 
value from WMAP. Lower panel: The total baryon fraction derived for the Vikhlinin et al. (2006) and Gastaldello et al. (2006) clusters if 
our best-fit stellar baryon relation is used to estimate a stellar baryon contribution for each of these systems. The weighted mean for the 
sample (dashed line) is Hfj/fim = 0.133 it 0.004. We observe no trend in baryon fraction with cluster mass. The error bars on the individual 
data points and the weighted mean include only statistical uncertainties. The unweighted mean for the combined Zhang et al. (2007) and 
Rasmussen & Ponman (2006) samples (dotted line) is included to provide a sense of the systematic uncertainties. 



3.5. Total Mass in the Intracluster Medium 

Determination of total cluster baryon fractions also re- 
quires measurement of the mass of hot gas in the intra- 
cluster medium. Because we lack these data for clusters 
in our sample, in §4 we fit the behavior of the stellar 
mass fraction with cluster mass and use this relation 
to derive total baryon fractions for clusters with pub- 
lished X-ray gas fractions. For the gas mass fractions 
we adopt the values of X-ray gas fraction inside of rsoo 1 
fg, and of enclosed mass, iWsooj from the tabulations by 
Vikhlinin et al. (2006) and Gastaldello et al. (2006). The 



latter adds only two new systems, but helps constrain the 
baryon fraction at low mass. We note that there are sev- 
eral systems for which the X-ray gas fraction was mea- 
sured by both groups. In these instances there is reason- 
ably good agreement, and hence we expect no significant 
relative systematic biases between the two samples. 

4. RESULTS AND DISCUSSION 

4.1. Baryon Mass Fraction versus Cluster Mass 

In this section, we obtain two principal results. First, 
summing the three baryonic components — the hot gas 
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of the intracluster medium, the stars of the BCG+ICL, 
and the stars of the other cluster galaxies gives a total 
baryon mass fraction, /{,, that is constant as a function 
of total mass (Figure 1). This constancy contradicts pre- 
vious studies and arises in large part from the inclusion 
of the ICL for systems with 51^500 < few xIO^^Mq. Sec- 
ond, the weighted mean value of the total baryon mass 
fraction is ft = 0.133.'* For reasons discussed below, we 
conclude that there is no compelling case for undetected 
baryons in the groups and clusters in our sample. 

We fit the behavior of the stellar and X-ray gas mass 
fractions with mass as power laws (top panel of Figure 1). 
Within rsoo, the stars result in a relationship log(/,) = 
(7.57 ± 0.08) - (0.64 ± 0.13) log(5M-5oo) and the plasma in 
log(/g) = (-3.87 ±0.04) + (0.20 ±0.05) log(«-5oo) for the 
sample of Viklilinin et al. (2006).''' fit the Viklilinin et al. 
(2006) sample alone The behavior of the two compo- 
nents is clearly inverted. The plasma component, which 
dominates at cluster masses > lO^^M©, has a shallower 
dependence on mass than the stellar component. We 
use our best fit relation for the stellar mass to convert 
the X-ray gas fractions from Vikhlinin et al. (2006) and 
Gastaldello et al. (2006) to total baryon fractions on a 
cluster- by-cluster basis. The resulting total baryon frac- 
tions, shown in the lower panel of Figure 1, are indepen- 
dent of cluster mass over the range of 5W500 from 6 x 10^^ 
to 10"'^^ Mq. If one considers only the statistical uncer- 
tainties associated with our total cluster baryon fraction 
and the WMAP measurement, then the two values are 
discrepant at the level of 3.2a-. 

There are three possible interpretations of this result in 
the context of a full baryon accounting. First, systematic 
uncertainties may be significantly larger than the random 
errors, and our measurement may therefore be consistent 
with the universal WMAP value. Because the gas mass 
fraction is much larger than the stellar mass fraction over 
the majority of the mass range we probe, we focus this 
discussion on potential errors in the gas mass fraction. 
Systematic errors, by their nature, are often difficult to 
calculate and can best be illuminated by comparing inde- 
pendent studies. Zhang et al. (2007) and Rasmussen & 
Ponman (2006) provide gas mass fractions for high and 
low mass systems, respectively. Zhang et al. (2007) find 
gas fractions that are 15% larger than those we adopted 
over the range probed by their sample. Rasmussen & 
Ponman (2006) also find gas fractions that are larger 
than those we adopted for comparably low mass systems. 
Together, these two samples cover the mass range of our 
adopted sample and result in an average /t, = 0.149, 
which is 2.1(7 discrepant with the WMAP value and 12% 
larger than the fi, we derive for the Vikhlinin et al. (2006) 
data set. We conclude that indeed the systematic uncer- 
tainties dominate, that the sense of the uncertainty works 
to diminish the discrepancy between measiirements and 
expectations, and that the case for undetected baryon 
components rests on removing these systematic uncer- 
tainties. Second, the baryon mass fraction with rsoo in 
clusters may be different from the universal value. Clus- 

Defining Y to be the ratio of the baryon fraction to the cosmic 
value from the WMAP third year results, this value for the total 
baryon fraction corresponds to Y = 0.76. 

Note that while this fit applies only to the Vikhlinin et al. 
(2006) data set, the Gastaldello et al. (2006) points in Figure 1 are 
consistent with this relation. 



tcr simulations routinely predict baryon fractions that 
are lower than the universal baryon fraction by ~10%. 
We have deliberately not adopted any of those as the 
target because there are still a minority of simulations 
that predict a baryon fraction larger than universal (e.g. 
Kravtsov et al. 2005). Third, the shortfall may be re- 
solved by yet undetected baryon components. The con- 
stancy of the measured baryon fraction with total system 
mass argues against a large undetected component, but 
one could contribute at a modest level and bypass current 
detection. In the end, all of these possibilities may con- 
tribute at some level to a final resolution of the baryon 
accounting, but at this point we are left to conclude that 
there is no compelling evidence for undetected baryons. 

We close this section by comparing our study to that 
of Lin et al. (2003), who obtained a dependence of the 
baryon fraction on cluster mass in conflict with our re- 
sult. The range of masses covered by the two studies is 
the same, and we both use rsoo as our fiducial radius, 
so the discrepancy is real. The most obvious differences 
are our inclusion of the ICL component and our different 
normalization of the gaseous component. We model the 
effect of including the ICL into their accounting by us- 
ing our relationship for the BCG+ICL mass vs. system 
mass and "correcting" their values of the baryon fraction 
on a cluster-by-cluster basis. ^ The relationship between 
baryon fraction and mass decreases from a result in 
the original Lin et al. (2003) study to a 1.8cr result in 
the corrected case. Therefore, the inclusion of the ICL 
explains part of the discrepancy between our result and 
that of Lin et al. (2003), although apparently not all of 
it. The lack of the ICL in the Lin et al. (2003) accounting 
was not an oversight — they attempted to model the ef- 
fect of the ICL in Lin & Mohr (2004), but simply lacked 
data of the sort we present here. While the exclusion 
of the ICL causes a relative underprediction of the con- 
tribution of stellar baryons for the lowest mass sytems, 
the Lin et al. adopted X-ray gas fractions, which are 
~ 25% larger than ours and hence consistent with the 
WMAP baryon mass fraction for the most massive sys- 
tems, cause a relative overprediction of baryons in the 
highest mass systems. These two effects lead to the ap- 
parent mass-dependent baryon fraction of Lin et al. As 
with the comparison among different studies discussed 
previously, we have an unresolved discrepancy in the gas 
mass fractions at the ~ 20% level that highlights the dif- 
ficulty in such an analysis and supports our contention 
that there is yet no compelling evidence for either addi- 
tional baryonic components or a trend in baryon fraction 
with system mass. 

4.2. Stellar Component of Clusters and Groups 

4.2.1. Behavior vs. 5W500 

We now consider the behavior of the stellar compo- 
nent in more detail. We compute the dependence of the 
optical luminosity on 51^500 for the cluster galaxies, the 
BCG+ICL, and the sum of these populations. There 
are multiple studies of the dependence of the total stel- 
lar luminosity of cluster galaxies on cluster mass. In 
the if— band, Lin et al. (2004) find that for the galaxies 
alone (excluding the BCG) L oc AT with 7 = 0.82 + 0.05 

^ This is a slight (~20%) overcorrection because it includes the 
BCG, which Lin et al. (2003) had already included. 
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1000 

log a [km sec"'] 

Fig. 2. — Calibration of velocity dispersion to X-ray derived 
masses within rgoQ. We plot the seven systems from Vikhlinin et al. 
(2006) for which we found velocity dispersions in the literature 
(Girardi et al. 1998; Wu et al. 1999). The line is our fit to the data, 
and hence our calibration for converting the velocity dispersions 
measured for our optical clusters into fWsoo's. 



(7 = 0.82 ± 0.04) within r5oo (7*200), and Ramella et al. 
(2004) find 7 = 0.74 ± 0.06 within raoo- A least-squares 
fit to our J— band data for the cluster galaxies yields 
7 = 0.71 ± 0.07 within rsoo, which is consistent with 
these results and most other published values (Girardi 
et al. 2000; Rines et al. 2004, but see Kochanek et al. 
(2003) for an exception). We conclude that the conflict 
between our result — that the baryon mass fraction is 
insensitive to system mass for groups and clusters — and 
those of past studies does not arise predominantly from 
our measurement of cluster galaxy properties. 

In comparison with the cluster galaxy luminosity 
versus cluster mass relation, the correlation between 
BCG+ICL luminosity and cluster mass is weaker, in- 
creasing by only about a factor of two in total luminosity 
between poor groups with a = 400 km s^^ and rich clus- 
ters with (7 = 1000 km s~^ (Figure 3). The existence of 
groups in which the fraction of stars in the BCG+ICL ex- 
ceeds that in galaxies out to rsoo (Figure 4) demonstrates 
that the rich cluster environment is not a requirement for 
the production of a significant intracluster stellar popu- 
lation. Instead, it appears that local density maxima 
in both the group and cluster environments are capa- 
ble of generating intracluster stars, presumably via tidal 
processes. This conclusion is valid whether one is inter- 
ested in the combination of the BCG and ICL, or in the 
ICL alone, because the latter contains the bulk (typi- 
cally > 80%) of the total light in the two components 
(Gonzalez et al. 2005). However, our results should not 
be interpreted to mean that all groups have a significant 
ICL component, because our original sample selection 
focused on systems with dominant BCGs. 

The observed weak dependence in the BCG+ICL is 
given by 7 = 0.39 ± 0.05. Physically, this result im- 
plies that production of the BCG+ICL — presumably 
from the stripping of stars from galaxies — is more effi- 
cient in group (or subcluster) environments than in clus- 
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Fig. 3. — The total absolute magnitude of the BCG+ICL com- 
ponent (circles) and of cluster galaxies (triangles) within rsoo as a 
function of cluster velocity dispersion. The filled circles denote sys- 
tems where the BCG+ICL contributes less than 40% of the total 
luminosity within this radius, while open-crossed circles correspond 
to systems with higher fractions of light in the BCG+ICL. The 
BCG+ICL luminosity depends less strongly on the velocity disper- 
sion than does the total galaxy luminosity. We also note that there 
is no systematic difference in the BCG+ICL absolute magnitudes 
for clusters with high and low BCG+ICL fractions. The BCG+ICL 
fraction differences are instead driven predominantly by variations 
in the total absolute magnitudes of the non-BCG cluster galaxy 
populations. 



tors. Again, we caution that when selecting clusters and 
groups for this sample, we specifically targeted systems 
with dominant BCGs. If there is a strong correlation 
between the presence of a dominant BCG (i.e., our sam- 
ple) and the total BCG+ICL luminosity (which is not 
necessarily the case because the ICL contributes much 
more light), then our selection criteria would artificially 
produce a smaller 7 than that actually found in nature. 

4.2.2. Total Cluster Luminosities and Mass-to-Light 

Ratios 

The observed weak dependence of the BCG+ICL lumi- 
nosity on cluster mass, if real, implies that 7 for the total 

stellar light within rsoo (isoo) niust be smaller than the 
value for the galaxy population alone. The best fit value 
for the behavior of L500 is 7 = 0.47 + 0.05, with a scatter 
of only 0.26 mag about the best fit. Recasting in terms 
of a mass-to-light ratio yields M500/L500 oc L^qq^^'"^^. 
For comparison, literature values that include only the 
cluster galaxy contribution typically have M/L (x L^ "^. 
Thus, if the slope we measure for the BCG+ICL com- 
ponent is not significantly biased by our exclusion of 
groups without a dominant BCG, then our data argue 
that galaxy cluster mass-to-light ratios are much more 
strongly dependent upon system mass than is commonly 
presumed. 

Conversely, if our selection against groups lacking a 
dominant BCG is the reason that the BCG+ICL and 
cluster galaxy luminosities have different observed clus- 
ter mass dependences, then there must exist a correlation 
between mass-to-light ratio and the presence of a dom- 
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Fig. 4. — The fraction of the total stellar luminosity, measured 
within rsoo, which is contained in the BCG+ICL component as 
a function of cluster velocity dispersion. The dotted vertical line 
corresponds to u = 750 km s~^, roughly the median velocity 
dispersion for an Abell richness class R > 1 cluster (ZabludoflF 
et al. 1990). 
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Fk;. 5. — The cumulative fraction of light contained in the 
BCG+ICL as a function of radius. Subdividing the sample 
into clusters with a > 750 km s^^ (open boxes) and a < 
750 km s^^ (filled triangles) demonstrates that the BCG+ICL 
fraction is systematically lower in the more massive clusters. The 
horizontal dashes included in the plots denote the rms variations 
for the cluster ensemble at each radius. 



inant BCG. Specifically, for a given velocity dispersion, 
groups with a dominant BCG have higher BCG+ICL lu- 
minosities and hence lower mass-to-light ratios. 

4.2.3. BCG+ICL Luminosity Fraction: Dependence on 
Velocity Dispersion and Radius 

We determine the fraction of total stellar luminosity 
contained in the BCG+ICL component as a function 
of velocity dispersion and radius. In Figure 4, we plot 
the BCG+ICL luminosity fraction, Lbcg+icl/ {Lgai + 
Lbcg+icl) a function of velocity dispersion within 
''500, whic-h is slightly larger than the physical region 
within which the BCG and ICL were modeled in Paper 
I. We find a trend of decreasing BCG+ICL luminosity 
fraction with increasing velocity dispersion. The same 
trend is observed if we use a fixed physical radius of 300 
kpc, but the scatter is higher. The presence of this trend 
is also independent of whether r5oo or a fixed fraction of 
r2oo is used for the abscissa. 

Within the clusters, the radial dependence of the 
BCG+ICL luminosity fraction is shown in Figure 5, 
where we plot the BCG+ICL luminosity fraction en- 
closed within different radii. Combining the full ensemble 
of clusters, the BCG+ICL luminosity fraction decreases 
smoothly from 65% within 0.1r2oo to 33% within r2oo- 
This radial decline, which indicates that the ICL is more 
concentrated than the galaxies, is also seen by Zibetti 
et al. (2005) and is predicted by the simulations of Mu- 
rante et al. (2004). 

Lastly, we use a threshold of ct = 750 km — 
roughly the median velocity dispersion of an Abell rich- 
ness class -R > 1 cluster (Zabludoff et al. 1990) to 
divide the data set into two subsamples (Figure 4). The 
clusters with higher velocity dispersions exhibit lower 
BCG+ICL luminosity fractions at all radii. To confirm 
that this result is not driven purely by the BCG, we also 
compute the corresponding luminosity fractions within 



circular annuli from to 0.5r2oo- We find that at all 
radii the fraction of the cluster luminosity contained in 
the ICL is lower for the more massive clusters. More- 
over, the rms scatter about the mean for the BCG+ICL 
luminosity fractions is only 10% for the high-cr sample 
compared with 20% for the low-cr sample within 0.5r2oo- 
Thus, for clusters and groups with a dominant BCG, 
the BCG+ICL luminosity fractions are smaller and more 
uniform in the more massive systems. 

4.3. Origin of the Baryon Trends 

We have identified two independent trends with mass 
that shed light on the environmental dependencies of 
physical processes that distribute baryons among the var- 
ious phases. First, we find that the fraction of stars in 
the ICL grows (and that in galaxies diminishes) with de- 
creasing system mass. This result is predicated on the 
assumption that we are not strongly biased in the type 
of system observed, as discussed above. Even so, we now 
know that in at least some low mass systems the ICL 
forms efficiently. As a result, we can rule out cluster- 
only processes, such as harassment (Moore et al. 1996, 
1998), as the dominant ICL formation mechanism, and 
find support for mechanisms that would be more effec- 
tive in lower mass systems, such as galaxy-galaxy inter- 
actions (Merritt 1984). Such tidal interactions, which 
could strip stars from galaxies and build up the ICL, are 
favored in lower mass groups, because the internal veloc- 
ity dispersions of group galaxies are similar to the veloc- 
ity dispersion of the group as a whole. If all low mass 
systems do have high ICL luminosity fractions, and some 
are the building-blocks of clusters, then a diluting mech- 
anism, such as the infall of individual galaxies into clus- 
ters, would be necessary to reduce the final ICL fraction 
of more massive systems. Of course, it may be that not 
all low mass groups have high ICL fractions like those in 
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our sample. The low ICL fractions of high mass clusters 
may therefore simply reflect their growth by the accre- 
tion of groups with lower ICL fractions than we observe. 
A related possibility is that the clusters accrete and mix 
progenitor groups before they have time to form substan- 
tial ICL via galaxy-galaxy interactions, whereas the ICL 
continues to grow in more isolated groups. Simulations 
of the formation of the ICL must reproduce these trends 
and will lead to a better understanding of the details of 
the formation mechanisms (Willman et al. 2004; Murante 
et al. 2004; Sommer-Larsen et al. 2005; Sommer-Larsen 
2006; Rudick et al. 2006). 

Second, we find that the baryon mass fraction across 
the range of systems for which we have both optical 
and X-ray data is roughly constant over the mass range 
6 X 10^3 < ^^^^ < 1 X 10^^ Mq. This result requires 
a fairly well-tuned balance between the X-ray gas mass 
fraction, which declines by a factor of two between sys- 
tems with !W5oo of 10^''^AfQ and 10^"* Mq, and the stel- 
lar mass fraction, which increases by a factor of four 
to compensate. The matching of the two trends sup- 
ports the suggestion (Bryan 2000; Lin et al. 2003) that 
the X-ray gas mass fraction behavior is due to the de- 
creased efficiency of turning gas into stars in the more 
massive environments. The reduction in star formation 
efEcicency in higher mass systems may arise from the 
decreased efficiency of tidal interactions among galaxies, 
the removal of the gaseous reservoirs of galaxies by inter- 
actions with the intracluster medium, or increased feed- 
back processess that quench star formation. As usual, 
interpreting such trends is complicated by the fact that 
today's lower mass systems are not necessarily similar to 
the lower mass antecedents of today's higher mass sys- 
tems. 

5. CONCLUSIONS 

We have conducted a study of 23 nearby galaxy clusters 
and groups with a dominant, early- type brightest cluster 
galaxy. For this sample we determine the luminosities of 
the different stellar constituents (including stars associ- 
ated with cluster galaxies and the intracluster stars) as 
a function of cluster velocity dispersion (cr) , cluster mass 
(iW'soo), and cluster-centric radius. Combining these re- 
sults with current measurements of the X-ray gas mass 
as a function of a, we calculate the dependence of the 
baryon budget on environment. We reach the following 
conclusions: 

1. The observed total baryon mass fraction within 
^■500 is constant for systems with iWsoo between 6 x lO-"^^ 
and 10^^ Mq. Although the mean baryon mass frac- 
tion is formally > 3a lower than the WMAP measure- 
ment, both the observational systematic uncertainties — 
which dominate the formal random errors — and the 
predicted deficit of baryons in clusters relative to the 
universal value are on the order of the observed baryon 
shortfall. Thus, we conclude that there is no compelling 
evidence for undetected baryons at this time. 

2. The inclusion of the combined light from the bright- 
est chister galaxy and the intracluster stars, which we 
refer to as the BCG-I-ICL component, becomes increas- 
ingly important in the baryon budget as the system 
mass decreases. For systems with 51/500= IO^'^Mq, the 
BCG-hlCL component is nearly 35% of the total stellar 
luminosity within r^oo, and the mass of the stellar com- 



ponent including the BCG+ICL and the cluster galaxies 
is comparable to that of the X-ray gas component. Con- 
sequently, the increase in X-ray gas mass fraction with 
increasing ^500 is directly reflected by a decrease in the 
total stellar mass fraction. The matching of the falling 
stellar mass and rising gas mass with total mass supports 
the suggestion made previously (Bryan 2000; Lin et al. 

2003) that the star formation efhciency decreases with 
increasing system mass. 

3. Considering only the stellar components, we ob- 
serve that the fraction of the total stellar luminosity in 
the BCG+ICL decreases with system velocity dispersion. 
The origin of this trend is that the luminosity of the 
BCG-I-ICL component changes more slowly with cluster 
mass than does the luminosity of cluster galaxies, in- 
creasing by only a factor of two between 400 kms~^ and 
1000 km s^^. The presence of large BCG+ICL luminos- 
ity fractions in the lower mass systems demonstrates that 
the rich cluster environment is not required for produc- 
tion of "intracluster" stars. 

The decrease in the importance of the BCG+ICL com- 
ponent, and in particular the dominant ICL component, 
with increasing cluster mass is consistent with more ef- 
ficient stripping of stars from galaxies - presumably via 
galaxy-galaxy interactions - in the lower mass systems. 
Although it remains unclear the degree to which our sam- 
ple selection, which favored systems with central, giant 
elliptical galaxies, impacts this result, it is noteworthy 
that Purccll et al. (2007) predict precisely the type of 
relation that we observe using a simple analytic model 
for BCG+ICL formation from satellite accretion and dis- 
ruption. A similar model is also able to reproduce the 
relative luminosities that we observe for the BCG and 
ICL (Conroy et al. 2007), implying that this basic physi- 
cal picture can largely explain two significant results from 
our work without invoking significant bias due to sample 
selection. 

4. The mean BCG+ICL luminosity fraction for our 
sample decreases monotonically from 65% within 0.1r2oo 
to 33% within r2oo. This trend is qualitatively similar to 
that seen in the stacked sample of Zibetti et al. (2005). 

5. Interior to r2oo, we find that the BCG+ICL optical 
luminosity fraction is inversely correlated with cluster 
velocity dispersion. Our mean value of 33% is some- 
what larger than one might expect based upon the work 
of Zibetti et al. (2005), who found a similar luminos- 
ity fraction of 32% within a smaller, 500 kpc radius. 
This difference likely reflects a difference between the 
systems included in the two studies. For clusters with 
a > 750 km s~^ , we obtain a mean value of 23% 
within r2oo, which is consistent with both the Zibetti 
et al. (2005) results and simulated BCG+ICL luminos- 
ity fractions within r2oo for rich clusters (Willman et al. 

2004) . 
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